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We argue that the a.c. conductivity a(oj) in the many-body localized phase is a power law of 
frequency io at low frequency: specifically, a(u>) ~ u> a with the exponent a approaching 1 at the 
phase transition to the thermal phase, and asymptoting to 2 deep in the localized phase. We 
identify two separate mechanisms giving rise to this power law: deep in the localized phase, the 
conductivity is dominated by rare resonant pairs of configurations; close to the transition, the 
dominant contributions are rare regions that are locally critical or in the thermal phase. We present 
numerical evidence supporting these claims, and discuss how these power laws can also be seen 
through polarization-decay measurements in ultracold atomic systems. 

PACS numbers: 


I. INTRODUCTION 


Isolated interacting quantum systems can undergo a 
dynamical phase transition—termed the many-body lo¬ 
calization (MBL) phase transition—between a “thermal” 
phase in which the system comes to thermal equilib¬ 
rium from generic initial conditions and a “localized” (or 
MBL) phase in which it does not [1-8]. Instead, an iso¬ 
lated system in the MBL phase is a “quantum memory”, 
retaining some local memory of its local initial condi¬ 
tions at arbitrarily late times [9-18]. The existence of 
the MBL phase can be proved with minimal assump¬ 
tions [20]; many of its properties are phenomenologi¬ 
cally understood [10, 11, 16], and can be explored us¬ 
ing the strong-randomness renormalization group [9, 21- 
23] . While the eigenstate properties of MBL systems are 
similar to those of noninteracting Anderson insulators, 
there are important differences in the dynamics, such as 
the logarithmic spreading of entanglement in the MBL 
phase [6, 9, 11, 18, 19, 24, 25]. 

In this work we show that the MBL phase also dif¬ 
fers sharply from noninteracting localized systems in its 
low-frequency response. We focus on a.c. conductivity 
in the MBL phase, for concreteness and to make con¬ 
tact with previous literature on solid-state systems (e.g., 
Mott’s law [26]); however, as we argue below, our dis¬ 
cussion directly extends to relaxation dynamics, which is 
more easily accessible in experimental studies using ul¬ 
tracold atoms [27], polar molecules [28], nitrogen-vacancy 
centers [29], and other forms of “synthetic” matter. We 
identify two physical mechanisms underlying the slow re¬ 
sponse: (a) the presence of resonant pairs of charge or 
spin configurations , connected by slow many-body rear¬ 
rangements; and (b) the presence within an MBL sys¬ 
tem of rare thermalizing regions, or “inclusions”, that 
act as local heat baths for their surroundings. These 
mechanisms are absent in noninteracting systems: thus, 
the differences in transport between single-particle and 


many-body localization can be traced to the much larger 
connectivity of the many-body Hilbert space. The two 
mechanisms we discuss involve local dissipative transport 
of the conserved densities (e.g., particle number), and are 
thus distinct from the “pure dephasing” processes that 
cause the slow growth of entanglement within the MBL 
phase [11, 25]. 

Our results for the a.c. conductivity are as follows. 
Whereas in noninteracting systems the a.c. conductivity 



FIG. 1: A.C. conductivity exponent a of a disordered one¬ 
dimensional system across the many-body localization tran¬ 
sition, showing four regimes: (i) the diffusive thermal phase; 
(ii) the subdiffusive thermal phase [31], which exists only in 
one dimension; and (iii, iv) the MBL phase. The MBL phase 
is divided into an “MBL-Griffiths” regime (iii) in which low- 
frequency response is dominated by rare critical or thermal 
regions, and an “MBL-Mott” regime (iv) in which it is dom¬ 
inated by pairs of resonant configurations. Our main pre¬ 
dictions are that the exponent a —► 1 [i.e., a ~ lo\ as the 
MBL transition is approached from the localized side, and 
that 1 < a < 2 throughout the MBL phase. These are 
consistent with numerical simulations of a nearest-neighbor 
random-field XXZ chain (shown in the plot). In the ther¬ 
mal phase, finite-size effects are strong, and a more careful 
analysis [31] is needed to extract the conductivity exponent. 
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a(co) ~ ui 2 log d+1 u) (Mott’s law [26] in d dimensions), in 
the interacting MBL phase at high temperature the con¬ 
ductivity goes as <j(u>) ~ w“, where a is an exponent that 
varies continuously throughout the MBL phase, ranging 
from a = 1 at the MBL transition to a —> 2 deep in 
the MBL phase. The exponent a has two regimes of be¬ 
havior, corresponding to the two mechanisms described 
above. Deep in the MBL phase, the conductivity is dom¬ 
inated, as in Mott’s law, by resonant transitions between 
localized configurations. It is enhanced relative to nonin¬ 
teracting localization because more such resonances are 
possible: in addition to single-particle hopping, a MBL 
state can undergo multiple-particle rearrangements. We 
term this regime the “MBL-Mott” regime, and argue that 
response in this regime is dominated by rare regions that 
are still localized but with anomalously large localization 
length. Close to the transition, the conductivity is dom¬ 
inated by rare thermalizing or critical regions and their 
surroundings; we call this the “MBL-Griffiths” regime. 
As the system approaches the transition from the insu¬ 
lating side, thermalizing inclusions proliferate; we show 
that this leads to the conductivity exponent a —> 1. 

In two or more dimensions, the exponent a = 0 
throughout the thermal phase (i.e., there is presumably 
a nonzero d.c. conductivity). In one dimension, however, 
a subdiffusive phase with a continuously varying conduc¬ 
tivity exponent 0 < a < 1 exists on the thermal side of 
the MBL transition [23, 30, 31]. Remarkably, therefore, 
the a.c. conductivity exponent a in one dimension is 
continuous and apparently smooth across the MBL tran¬ 
sition, approaching the critical behavior a(co) ~ w from 
both sides (Fig. 1). 

This paper is structured as follows. In Sec. II we 
list our assumptions. In Secs. Ill and IV we introduce 
the MBL-Mott and MBL-Griffiths phases respectively; 
then in Sec. V we discuss the transition between these 
phases. In Sec. VI we describe the numerical methods 
used to compare our theoretical predictions with data 
on random-field Heisenberg spin chains (details of the 
numerical analysis are given in Appendix B). Sec. VII 
connects the a.c. conductivity to the relaxation dynam¬ 
ics measured in ultracold atomic experiments. Finally, 
Sec. VIII summarizes our results and comments on their 
scope. 


II. ASSUMPTIONS 

We work with a generic disordered lattice Hamiltonian 
having a conserved density (e.g., a particle number, or 
a particular projection of spin). The current associated 
with this charge is denoted j. The a.c. conductivity 
tensor a in the T —y oo limit is then given by the Kubo 
formula: 

Ta a0 (uj) = — E( m\j a \n)(n\jp\m)5{u - u) mn ) . (1) 


Here T is the temperature, N is the number of sites, Z 
is the partition function which in this infinite T limit is 
the dimension of the many-body state space; the indices 
in, n run over all Z many-body eigenstates; and the cur¬ 
rent j a is the sum over local currents, viz. j a = ]>//■ ji. a - 
We shall only be concerned with the diagonal elements 
cr aa , so henceforth we shall drop the index a. Our argu¬ 
ments should also apply to the frequency-dependence of 
the a.c. thermal conductivity, e.g., in systems where the 
only conserved quantity is the energy. 

When we consider the MBL phase, we specialize here 
to the case where all many-body eigenstates are localized, 
so we can discuss in terms of the localized conserved op¬ 
erators. However, the results we obtain should also apply 
to the MBL phase in systems with a many-body mobility 
edge. In the latter case, when we discuss ‘rare regions’ 
they are not only rare local disorder configurations in the 
system’s Hamiltonian, but also rare local configurations 
of the state that put it locally closer to, at, or across the 
mobility edge. 

Our considerations here are at the level of linear re¬ 
sponse theory: i.e., we assume throughout that the drive 
is sufficiently weak and is present for a sufficiently short 
time that linear response applies. It has recently been 
shown [19] by one of the present authors that local¬ 
ized systems subject to a fixed-amplitude drive display a 
highly non-local response at low enough frequencies. Fur¬ 
ther, an MBL system subject to a finite-frequency drive 
for a long enough duration will eventually leave the linear 
response regime and enter instead a Floquet MBL steady 
state or even thermalize due to the ac drive [8, 19, 32]. 
In the present work we are not concerned with these 
regimes. 

III. MBL-MOTT REGIME 
A. Many-body “Mott” conductivity 

We begin by considering the generic behavior deep in 
the MBL phase. As discussed earlier, we specialize to the 
regime where all eigenstates of the system are localized. 
In this regime, the system Hamiltonian admits a repre¬ 
sentation in terms of effective spin-1/2 degrees of freedom 
labeled r^, which are frequently referred to as local inte¬ 
grals of motion or “1-bits” [10, 11, 16]: in terms of these, 
II = XZ, hiT? + J2i,j Jij T i T j + .... Eigenstates of II are 
also simultaneously eigenstates of all the r^. These ef¬ 
fective r spins are related to the microscopic degrees of 
freedom (which need not be spin-1/2) by a unitary trans¬ 
formation that is local up to exponentially small tails. In 
terms of the effective r spins, the current operator can 
be expressed as 

J = E <K Q + E + ... , (2) 

a.,k 

where a = x,y, z\ r“ is the appropriate Pauli matrix; 
and k , l run over effective spins. The coefficients 
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for n > 2 fall off exponentially with the distance between 
the farthest effective spins in that term. Stability of the 
MBL phase further requires them to fall off exponentially 
with the order n [33] . Note that for a single-particle (non¬ 
interacting) Anderson insulator the r“ operators are the 
creation, annihilation and number operators of the lo¬ 
calized single-particle states, so the coefficients K are 
zero for n > 2 , i.e., the current operator only contains 
single-particle hops and no multiple-particle rearrange¬ 
ments. 

We now briefly review Mott’s argument [26] for the 
a.c. conductivity in noninteracting localized systems at 
temperature T > ui. For this noninteracting case, the 
transitions contributing to <t(oj) at low frequency involve 
rare pairs of resonant sites that hybridize to form pairs 
of nearly-degenerate eigenstates (i.e., symmetric and an¬ 
tisymmetric combinations of the wavefunctions centered 
at the two sites) with small energy splitting to. Short- 
distance resonances, while common, typically have large 
splittings because of local level repulsion; these local pro¬ 
cesses give only a subdominant contribution to the con¬ 
ductivity in the low-frequency limit [34]. To find the 
resonant pairs of sites with energy splitting ui that dom¬ 
inate in ct(w), one has to go a distance r u determined 
by the condition W exp (—r w ) = w, where W is a micro¬ 
scopic bandwidth. (We are assuming the single-particle 
localization length is of order one lattice spacing and do 
not include factors of it.) The number of such pairs is 
~ r'fifi 1 . The typical current matrix element between 
such pairs of eigenstates is wr u ~ u>\og(W/uj) (because 
they involve moving a unit charge a distance r u at a rate 
w). Putting these pieces together, we recover the Mott 
result: 

er(w) ~ w 2 \og d+1 {W/u >). (3) 

Note that the contribution from more distant pairs is 
weaker, because the current matrix element falls off as 
exp(— R) whereas the phase space only grows as 

This argument is fundamentally altered by many-body 
processes for the interacting MBL phase. Here the con¬ 
ductivity includes not only hopping resonances between 
pairs of sites but also many-body resonances between 
pairs of configurations ; hence the “phase space” factor is 
strongly enhanced. This argument is fundamentally al¬ 
tered by many-body processes for the interacting MBL 
phase. Here the conductivity includes not only hopping 
resonances between pairs of sites but also many-body 
resonances between pairs of configurations. Hence the 
“phase space” factor is strongly enhanced. We now ar¬ 
gue that this enhanced phase space factor grows expo¬ 
nentially in the number of effective spins flipped. 

The many-body resonances that dominate the low- 
frequency dynamics flip n effective spins, with those spins 
typically having random spacings of order the localiza¬ 
tion length or less so that they do interact with each 
other. For d > 1 this set of spins will in general have 
a fractal geometry. Let 7 collectively denote all the rel¬ 
evant parameters (shape, typical interparticle spacing, 


etc.) specifying the ensemble of possibly resonant ‘clus¬ 
ters’ of flipped spins. (Given a cluster, in other words, 
one can characterize it through its parameters 7 ; differ¬ 
ent resonant n-spin clusters with the same 7 will have 
the same hybridization strength.) 

The typical current matrix element for a reso¬ 
nance with parameters 7 that flips n spins is ~ 
Wexp(—n/C( 7 )).Here (( 7 ) is a dimensionless quantity 
that depends on 7 and varies continuously in the MBL 
phase; £( 7 ) remains finite at the MBL transition and de¬ 
creases as one moves deeper into the localized phase. £( 7 ) 
is larger for resonances having more closely spaced and 
thus more strongly interacting spins.Let us fix these pa¬ 
rameters 7 . Then, analogous to the single-particle case, 
the frequency to picks out an “optimal” n such that 

W exp (—n/£(y)) = cu => n = (( 7 ) log (W/ 10 ). (4) 

The number of possible resonances (in the ensemble pa¬ 
rameterized by 7 ) that flip n effective spins in the im¬ 
mediate vicinity of one particular real-space location is 
exponential in n, while the frequency bandwidth of such 
rearrangements is linear in n. Thus, to leading order the 
density of states of possible resonances at order n grows 
exponentially with n. Specifically, it grows as ~ e s ^ n : 
where s(y) is the configurational entropy per flipped spin 
of the possibly resonant clusters in the ensemble 7 . This 
is the entropy of all the possible choices of the n spins 
flipped by the resonance. Using this and Eq. (4), the 
density of states of resonant configurations from ensem¬ 
ble 7 at frequency ui grows as a power law, where 

<j) = s( 7 )C( 7 ), i n contrast with the logarithmic growth 
in the noninteracting case. The dominant resonances at 
low frequency flip many spins and have their properties 
7 chosen so that the product s( 7 )£( 7 ) is maximized [35]. 

We now assume that we have maximized this product 
0, and complete our estimate of the MBL-“Mott” con¬ 
ductivity. The current matrix elements remain ~ w, up 
to logarithmic factors. Putting this together with the 
phase space factor w -< ^, the conductivity goes as 

er(w) ~ U] 2 ~^ (5) 

at low frequency. 

For the MBL phase to be stable, we need that a typi¬ 
cal eigenstate is, at a typical real-space location, not in¬ 
volved in many resonances. From the discussion above, 
the typical accessible phase space for final states with a 
matrix element of oj goes as ~ thus the typical level 
spacing for these goes as ~ oA In order that long-range 
resonances remain rare and do not destabilize the MBL 
phase, the matrix element must vanish faster than the 
typical level spacing in the long-distance, small to limit. 
Thus 0 < <fi < 1 (and thus 2 > a > 1) in the MBL phase, 
with (f> increasing (thus a decreasing) as the phase tran¬ 
sition to the thermal state is approached. 
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B. Rare-region Mott resonances 

In the above discussion, we argued that the low- 
frequency conductivity in the MBL phase is dominated 
by rare many-spin resonances, and goes as oj 2_s< ”, where 
s and ( are properties of the MBL phase (optimized over 
families of resonances parameterized by 7 ). However, in 
a disordered system, ( is itself a random variable, so there 
will be atypical clusters in which (for example) the ran¬ 
dom fields are small and therefore the system is locally 
closer to the delocalization transition. (We focus on ( 
but the same argument can be applied to any other pa¬ 
rameter.) In such segments, ( will take a local value Cioc. 
that deviates from its typical value (, and the matrix el¬ 
ement for resonances involving n spins will be atypically 
large. 

These rare local “regions” occur with a probability 
~ exp[—r/(Ci oc .)], where /(Cioc.) is a nonnegative “rate 
function” [36] that vanishes quadratically at Cioc. = C- 
By the above arguments the contribution of such a rare 
local resonance to the ac conductivity will be ~ w 2 , 
while the number of spins flipped by the rare resonance 
is n ~ Cioc. log(W/w) (which sets a minimum “size” for 
the rare region). Therefore the density of such rare reso¬ 
nances will be ~ w (Ci°c./(Cioc.)-Cioc.s)^ B ecause the initial 
‘gain’ in conductivity by going to these rare local reso¬ 
nances is linear in (Cioc. — C) while the probability ‘cost’ is 
only quadratic in the deviations from typical, the generic 
situation in the MBL phase is that the low-w conductiv¬ 
ity is dominated by rare many-body resonances in rare 
regions that are locally atypically close to the delocaliza¬ 
tion transition (i.e., have an atypically large Cioc.)- 

When the system is deep in the MBL phase, the domi¬ 
nant contributions to the low frequency conductivity are 
from resonant clusters in regions that are themselves in 
the localized phase; we call this regime the “MBL-Mott” 
regime. In the low frequency limit in this regime, each 
resonant cluster is large compared to its local value of 
the localization length. As the transition to the thermal 
phase is approached, at some point before reaching the 
transition these dominant rare clusters become instead 
locally critical or thermal quantum Griffiths regions. We 
now turn to such rare-region Griffiths effects, and show 
that they give rise to a conductivity exponent that ap¬ 
proaches a = 1 at the critical point. 


IV. MBL-GRIFFITHS REGIME 

We shall eventually be concerned with both thermal 
and critical rare regions, but to set up our discussion we 
begin by considering an inclusion that is locally deep in 
the thermal phase, embedded in a typical insulating en¬ 
vironment. This thermal inclusion is of volume V, has a 
many-body level spacing A that decreases exponentially 
with V and a transport time (i.e., Thouless time) txh 
that increases polynomially with V (specifically, as V 2 ^ d 
for a compact internally diffusive inclusion and with a 



FIG. 2: Left: schematic of the structure of a thermal inclusion 
in the MBL phase, showing its core (the microscopically rare 
thermal region), the periphery (the typical surroundings that 
get strongly entangled with the rare region), and the typical 
MBL surroundings where the core does not flip the local effec¬ 
tive spins. Low-frequency transport occurs through transport 
between the periphery and the core. Right: numerically com¬ 
puted conductivity exponent a of an XXZ chain of size L = 12 
with integrability-breaking next-to-nearest neighbor exchange 
coupling that contains a thermal or critical inclusion (i.e., a 
region that is locally thermal or critical) of four sites. IFtyp 
is the typical value of disorder, and three different disorder 
values for the inclusion Wi nc are shown. When the typical 
system is in the MBL phase (Ht yp > 8), the conductivity ex¬ 
ponent saturates to near one, consistent with the discussion 
in the main text. 


larger power for a fractal or critical inclusion). In gen¬ 
eral, < 1 for large thermal inclusions. (We ex¬ 

pect this also to be true for critical inclusions in d > 1 .) 
Moreover, each inclusion thermalizes its immediate insu¬ 
lating surroundings. Thus the inclusion consists of two 
parts: first, the ‘core’ of the inclusion which is the rare 
region that is locally thermalizing (or critical), and sec¬ 
ond, the typical insulating region surrounding this core, 
which gets strongly entangled with the core in the many- 
body eigenstates (we call this the ‘periphery’). A more 
thorough discussion of these inclusions is presented in 
Appendix A. 

Now we consider the a.c. response of this inclusion, 
probed at a frequency ui ~ A <C 1/trh- Because the core 
relaxes rapidly compared with co, it essentially adiabati- 
cally follows the applied electric field, and its response is 
reactive rather than dissipative [38, 39]. Specifically, to 
leading order, the core response goes as cr(w) ~ uj 2 tTh, 
which is subleading at low frequencies to the many-body 
Mott contribution. Thus transport within thermal inclu¬ 
sion cores does not dominate the low-frequency conduc¬ 
tivity. 

However, the periphery of an inclusion with core level 
spacing A does contribute strongly to its conductivity at 
frequencies down to A, as we now argue. This periphery 
consists of typical MBL regions that experience the core 
as a finite bath [40, 41] to which they are coupled with 
matrix elements that fall off as ~ exp(— R/() (where R 
is the distance from the core and C is a decay length). 
We can estimate the decay rate of a peripheral spin, us¬ 
ing the Golden Rule, as 7 (R) ~ W exp(—2 R/(). So long 
as 7 (R) A, the Golden Rule is valid on these time 
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scales and the core does indeed act as a ‘bath’ for these 
spins. Far from the inclusion core, however, 7 (R) <C A; 
a spin at this distance resolves the discreteness of core 
levels and does not decay into them. The overall pic¬ 
ture is as follows (Fig. 2): The core (with level spacing 
A) is surrounded by ‘shells’ of continuously decreasing 
7 (R), with the outermost ‘active’ shell having a decay 
rate j(R) — A. Beyond this distance the system is in¬ 
sensitive to the presence of the thermal core and remains 
fully localized. 

This picture thus gives the behavior of the a.c. conduc¬ 
tivity in the presence of a single such thermal inclusion. 
When one probes the system at a frequency u> > A, the 
conductivity is dominated by the shell at radius R u such 
that 7 (Ru) = u). Shells closer to the core relax faster, and 
their response to a probe oscillating at ui is mainly reac¬ 
tive; meanwhile, shells that are farther do not respond 
at all at u>. The conductivity of the dominant shell is 
proportional to its Golden-Rule decay rate, so this shell 
gives a ~ 'y(Ru) ~ ui (up to logw factors due to size, 
dipole matrix element, etc.). The conductivity due to a 
single inclusion thus turns on at and above a frequency 
A and has the behavior <7 ~ u> at intermediate frequen¬ 
cies. This reasoning extends directly to any inclusion 
whose core has an internal relaxation rate greater than 
its many-body level spacing, and is supported by numer¬ 
ical simulations (Fig. 2b) on inhomogeneous systems, in 
which thermal or critical inclusions are put in by hand: 
the conductivity in the presence of an inclusion goes as 
a ~ ui whenever the core is thermal or critical, as ex¬ 
pected. 

We now use the single-inclusion result to study the 
rare-region contribution to the a.c. conductivity of a 
generic MBL system, which contains some density of in¬ 
clusions at all scales. The cores of these inclusions must 
be thermal, but there does not appear to be any con¬ 
straint on how thermal they are, so the most common 
such cores of a given size will be cores that are locally 
arbitrarily close to the critical point. The inclusion cores 
that dominate the conductivity at low w are thus rare 
locally critical regions with level spacing u) and conse¬ 
quently of volume ~ Cc log(W/w); the probability of such 
cores is therefore ~ pCdog{w/u) ^ w j iere p j s (heuris- 
tically) the probability that a unit-volume region is lo¬ 
cally critical. (One can define p more precisely as follows: 
the density of critical inclusions of volume V decreases 
exponentially with V, as p v .) In the MBL phase p < 1, 
and p approaches one at the transition; thus the Grif¬ 
fiths exponent g , which is positive, approaches zero as 
the transition is approached. Since each such inclusion 
contributes ~ w (up to logarithmic corrections) to the 
conductivity, the resulting conductivity of the Griffiths 
insulator goes as cr ~ w 1+s , where the Griffiths exponent 
g goes to zero at the critical point and rises smoothly in 
the MBL phase. 

We briefly comment on how these Griffiths arguments 
connect with those in the thermal phase [23, 31]. In the 
MBL phase, as discussed above, thermal inclusions of 


large volume V become exponentially rare in V. On the 
thermal side, instead, it is localized inclusions that be¬ 
come exponentially rare at large scales. In one dimen¬ 
sion, these rare localized inclusions act as transport bot¬ 
tlenecks, leading to a subdiffusive Griffiths phase [31]. 
In higher dimensions, however, rare insulating regions in 
the thermal phase cannot block transport, and the d.c. 
conductivity is nonzero. 

V. TRANSITION BETWEEN MBL-MOTT AND 
MBL-GRIFFITHS REGIMES 

The overall behavior of the low-frequency conductiv¬ 
ity exponent a is shown in Fig. 1: everywhere in the 
MBL phase 1 < a < 2. Near the MBL phase transi¬ 
tion, critical and thermal inclusions proliferate and the 
dominant mechanism is “Griffiths”; far from the MBL 
transition, such inclusions are too rare, and the domi¬ 
nant contributions to conductivity are instead from rare 
many-body Mott resonances within locally insulating re¬ 
gions (which are still less insulating than the typical re¬ 
gion). These regimes transition into each other as fol¬ 
lows: Within the “MBL-Mott” phase the dominant re¬ 
gions are still locally insulating; as one moves towards 
the transition, these dominant regions become less in¬ 
sulating. Eventually, before the MBL phase transition, 
the dominant regions become critical, and the system en¬ 
ters the “MBL-Griffiths” phase. Throughout the MBL- 
Griffitlrs phase, the dominant rare regions remain critical, 
and only their prevalence changes: as the critical point 
is approached, these rare critical regions become more 
common, and eventually proliferate. Thus the physics 
underlying the evolution of the conductivity exponent is 
qualitatively different in the two regimes, and we expect 
that the exponent is nonanalytic (though perhaps quite 
smooth) at the “Mott-Griffiths” transition between these 
two regimes. The location of the Mott-Griffiths transi¬ 
tion line in Fig. 1 is schematic—determining the location 
of this transition within the MBL phase remains an in¬ 
teresting direction for future work. 

VI. NUMERICAL SIMULATIONS 

We have checked these expectations against simula¬ 
tions of the conductivity in the random-held XXZ chain, 
governed by the Hamiltonian 

h = \ J2{sts~+h.c.)+j z i+E > ( 6 ) 

(ij) * * 

where hi is a local quenched random held drawn uni¬ 
formly from [— W, W], J is the spin exchange energy 
scale, and J z the spin-spin coupling strength. We mea¬ 
sure energies in units of J; in all presented data, we also 
take J z = J = h = 1. We diagonalize the full Hamilto¬ 
nian to calculate the conductivity <r(u) at inhnite tern- 
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FIG. 3: Main panel: conductivity as a function of frequency in 
the random field XXZ chain for system size L — 14 and disor¬ 
der bandwidth W = 8, showing multiple decades of power-law 
behavior. Inset: behavior of the conductivity exponent deep 
in the MBL phase is consistent with the many-body Mott 
prediction (see main text) that 2 — a ~ £, combined with the 
perturbative estimate £ ~ l/log(IV/J). 


perature using Eq. (1) by binning it on a logarithmi¬ 
cally spaced frequency grid which typically ranges from 
(jj = 10 -6 to 2. In Figs. 1 and 3 we have only nearest- 
neighbor exchange, while we also have second-neighbor 
exchange of strength ,/' = 1 in Fig. 2. The rationale 
for the latter choice is that, while next-nearest neighbor 
exchange interactions exacerbate finite-size effects, they 
also break integrability and make the disorder-free delo¬ 
calized phase behave thermally, even for relatively small 
system sizes. Thus they are essential for correctly simu¬ 
lating small thermal inclusions with W[ nc = 0. 

The conductivity exponent a, shown in Fig. 1, lias 
been extracted from power-law fits to the low frequency 
response that hold over multiple decades, see Fig. 3 main 
panel for the example of W = 8 which yields a ~ 1.5. 
The numerical results further confirm that the conduc¬ 
tivity exponent a = 1 at the MBL transition and that 
it asymptotes to a ~ 2 in the strong disorder limit, see 
inset of Fig. 3. A more detailed discussion of some nu¬ 
merical issues and a comparison between noninteracting 
and many-body insulators is given in Appendix B. 


VII. EXPERIMENTAL ASPECTS 

The predictions in our work concern the finite-time 
dynamical properties of MBL systems; thus, they are ro¬ 
bust against weak coupling to an external bath, which 
is present in all physical systems [40]; so long as w > T, 
where T is the bath-induced linewidth, the bath will not 
change these conductivity power laws. Therefore our pre¬ 
dictions for conductivity can be tested experimentally, 
both in solid-state systems [42] and in ultracold atomic 
systems [27]. In electronic systems, a.c. conductivity is 
straightforward to measure, but the long-range nature of 


the Coulomb interaction will modify several of our con¬ 
clusions. 

In principle one can also measure a.c. conductivity 
in ultracold atomic systems such as optical lattices by 
applying a periodically modulated tilt to the entire lat¬ 
tice [43]. However, in the current optical-lattice MBL 
experiments [27] it is more convenient to study relax¬ 
ation in the time domain; we now show how our results 
generalize to such experiments. (Note: Ref [27] used a 
quasiperiodic potential, while we are discussing the case 
of a random potential.) In general, these experiments 
involve creating a particular nonequilibrium density con¬ 
figuration and measuring the evolution of its “contrast” 
(i.e., the overlap between final and initial density devia¬ 
tions from thermal equilibrium). In the MBL phase, this 
contrast (which we denote C(t)) approaches a nonzero 
saturated value (Ax,; we argue that it does so at long 
time t as 

C(t) - Coo - * 1_a , (7) 

where a > 1 is the a.c. conductivity exponent discussed 
here and plotted in Fig. 1. This result holds in both 
the MBL-Mott and MBL-Griffiths regimes. In the Mott 
regime, at time t, resonant pairs with splitting < 1/t 
are still in their initial state and retain their initial den¬ 
sity deviation, whereas faster pairs oscillate and thus 
have “forgotten” their initial density deviation. Counting 
all Mott pairs with splitting < ui, using the arguments 
above, we find that these go as = w“ _1 , which gives 
Eq. (7). Likewise, in the Griffiths regime, the contribu¬ 
tion at time t is due to the peripheral spins of inclusions 
with core level spacing A < 1/t. The density of such in¬ 
clusions is ~ A“ _1 , which once again yields the result (7). 
Note that this decay becomes very slow as the transition 
is approached: the exponent (l — a) goes to zero at the 
transition. These arguments apply to the long time be¬ 
havior when C (t) is near Coo; the earlier time regime near 
the critical point when Coo is small should be governed 
by the dynamical critical behavior. Preliminary numeri¬ 
cal simulations on XXZ chains suggest that, deep in the 
MBL phase, this rare-region contribution might be dif¬ 
ficult to detect in experiment, as its amplitude is small 
compared with steady-state fluctuations of the contrast 
in finite systems. 


VIII. DISCUSSION 

In this paper we have argued that the low frequency 
a.c. conductivity in the MBL phase goes as a(u>) ~ w a , 
with 1 < a < 2 throughout the phase, and a —> 1 as 
the delocalization transition is approached. Deep in the 
MBL phase, the dominant processes involve transitions 
between rare configurations, in rare regions that are lo¬ 
calized but have an anomalously large localization length. 
Near the transition, the dominant rare regions are locally 
thermal or critical instead. The power-laws we expect on 
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general grounds are consistent with those seen in numer¬ 
ical results for the random-field XXZ model. We em¬ 
phasize that the power laws we find in the optical con¬ 
ductivity are not related to those predicted for electron 
glasses [44] : we are considering high-temperature behav¬ 
ior (i.e., w«T) in models with short-range interactions, 
whereas those works consider low-temperature behavior 
(i.e., to T) in models with Coulomb interactions. 

We conclude with some comments on the scope of our 
results. As already discussed above, our analysis of a.c. 
transport directly extends to relaxation dynamics. More¬ 
over, our results here should also describe, e.g., thermal 
transport, in systems where the only conserved quantity 
is the energy. However, our analysis of the conductivity 
relies on the fact that the conductivity is related to the 
spectral function of a current (i.e., a quantity associated 
with a globally conserved charge) and does not extend 
to generic spectral functions, such as those probed using 
optical lattice modulation spectroscopy [43]. 

Our analysis also depends on both the disorder corre¬ 
lations and the interactions being short-range. Specifi¬ 
cally, we assume that the effective interactions that me¬ 
diate many-body resonances fall off exponentially with 
distance and with the number of effective spins involved. 
Thus our conclusions are modified in an essential way 
when the interactions instead fall off as a power law of 
distance; this case will be treated elsewhere. It is not 
presently clear whether or not the stretched-exponential 
effective interactions that occur at putative critical points 
within the MBL phase [9, 14, 21, 22] substantially modify 
the above story. Also, our analysis of near-transition be¬ 
havior assumes that the delocalized phase is thermal, and 
thus may not apply to hypothesized transitions between 
an MBL phase and a nonthermal delocalized phase [45- 
47]. 
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Appendix A: Structure of a thermal inclusion 

A thermal inclusion core (i.e., a large region with rare 
microscopic parameters) in the MBL phase acts as a lo¬ 
cal, discrete “bath” for the peripheral insulating mate¬ 
rial around it. Thus, in the many-body eigenstates it is 
strongly entangled with the nearby (“peripheral”) typi¬ 
cal regions. Because the inclusion is finite, sufficiently far 
from it the MBL phase with area-law eigenstate entan¬ 
glement re-establishes itself. In this Appendix we discuss 
how this crossover takes place. 

A naive estimate (which will turn out to be largely 
correct) is as follows: If one ignores the discreteness of 
the bath levels, a Golden Rule estimate [40] suggests that 
the decay rate of a typical degree of freedom (which for 
convenience we shall call a spin) a distance R from the 
inclusion core is ~ exp(—2 R/Q where the decay length ( 
remains finite as the MBL transition is approached. This 
rate must be compared with the many-body level spacing 
A ~ e~ sV of the inclusion core; when this putative decay 
rate is smaller than the many-body level spacing, the 
inclusion is actually unable to act as a bath [40] and no 
decay takes place. Thus, an inclusion core of volume V 
is surrounded by a thermal periphery of linear size ~ V. 
The characteristic relaxation rate decays from its value 
at the center of the core (which is of order the bandwidth 
for a strongly thermal inclusion) to its value at the edge 
of the periphery, which is ~ A. 

However, this argument is evidently incomplete. When 
the inclusion core thermalizes a spin, on sufficiently long 
time scales this additional spin is also “thermal”, and 
thus naively might be thought able to act as a bath for 
other, more distant, spins. If one iterates this reasoning, 
however, one arrives at an obviously incorrect result: an 
inclusion core of linear size L thermalizes a region of lin¬ 
ear size ~ L d around it; and the combined level spacing 

of this full thermalized region is now ~ e~ sL (s being 
the thermal entropy per spin), which naively allows it to 
thermalize yet further regions, and so on, until the entire 
MBL system is thermalized. To avoid this conclusion, 
one must understand why these peripheral regions that 
are “thermalized” by the core cannot act as a bath for 
more distant insulating regions. 

One can see this as follows: Let us first remove all cou¬ 
plings that cross the boundary between the thermal in¬ 
clusion core and the periphery. Then the Hamiltonian of 
the now MBL peripheral region can be written in terms 
of 1-bits, in terms of which it takes the fully diagonal 
form H = JY hiT? + JijT^T^ + ...; the thermal inclusion 
core is, of course, described by a generic thermalizing 
Hamiltonian. Now we reinstate the boundary couplings; 
these are local in terms of the physical spins, and thus 
generally consist of a strictly local physical operator O 
on the MBL side of the boundary, coupled to an oper¬ 
ator on the thermal side. The operator O involves 1-bit 
flips at all distances, but contributions from distant 1-bits 
are exponentially suppressed. Because the intrinsic 1-bit 


Hamiltonian is purely diagonal, an 1-bit at a distance particular, the nearer 1-bits do not act as a bath, and the 
l £ from the boundary can only thermalize through ability of an inclusion to tlrermalize its surroundings is 
its exponentially weak contribution to the operator O; in determined by the size of its core. 


Appendix B: Details of conductivity numerics 

In this section, we discuss some of the subtler issues involved in numerically extracting the conductivity exponents. 
We begin with a discussion of finite-size effects and boundary conditions. We then compare our many-body a.c. 
conductivity numerics with a study of single-particle (i.e., noninteracting) insulators at similar system sizes. We find 
that the size- and disorder-dependence of the many-body conductivity is consistent with theoretical expectations, and 
qualitatively different from that of the single-particle conductivity. Finally, we present data elucidating the nature of 
transitions contributing to the low-frequency conductivity. 


a. Boundary conditions and finite-size effects 

In numerics on finite systems, the conductivity exponents discussed here only occur at intermediate frequencies, 
wl C w < J, where wl ~ exp {—L/C,) in the MBL phase [or exp(—L/£) in the single particle case] is a size-dependent 
low-frequency cutoff, ut sets the scale for level repulsion between states or configurations that differ on length scales 
on the order of the system size. The behavior below this frequency scale depends on whether the boundary conditions 
are open or periodic. In the case of open boundary conditions, the conductivity at the lowest frequencies goes as 
cr(uj) ~ w 3 . As in the Mott argument, two factors of ui are due to the current matrix element, which is constrained by 
the boundary conditions to vanish as to at low frequencies. The third factor is due to level repulsion in the Gaussian 
Orthogonal Ensemble [39], and captures the phase space of pairs of states with these very small energy differences. 
For periodic boundary conditions, on the other hand, the conductivity in the finite-size-dominated regime scales as 
u(tu) ~ ui. In that case, there is still a factor of u> from level repulsion, but the current matrix elements do not vanish 
in the limit of low frequencies, because these currents “wrap around” the system and hence do not build up large 
charge imbalance even when they are at very low frequency. (The distinction between the two kinds of boundary 
condition can be intuitively understood by contrasting the behavior of (a) a finite metallic grain embedded in an 
insulator and subject to spatially uniform a.c. electric field, and (b) a conducting ring with an a.c. magnetic flux 
through it. The response of the former becomes essentially dielectric in the limit of low frequencies, whereas that of 
the latter remains dissipative.) 

This finite-size dominated regime is clearly seen in numerical simulations on the many-body system in the thermal 
regime [Fig. 4(a)]. Deep in the localized regime, cr(w) is insensitive to boundary conditions in the frequency range 
we can access [Fig. 4(c)]; the finite size effects presumably appear only at even lower frequency u>l ~ exp {—L/C,). In 
the near-critical regime (IF = 4), where a is near one, the finite-system behaviors for open and periodic boundary 
conditions are qualitatively different [Fig. 4(b)]: For periodic boundary conditions, the low frequency regime due to 
the finite size effect is expected to also have a = 1. No finite size effect is apparent at IF = 4 for periodic boundary 
conditions, which is possibly a consequence of the finite size regime having essentially the same scaling as the “bulk” 
regime. On the other hand, for open boundary conditions the low frequency finite size regime will have an effective 
exponent a e s = 3. The beginnings of the crossover in to this regime are apparent in Fig. 4(b), and we can see the 
results are converging with increasing L towards the periodic boundary condition results. 

The analogous results for noninteracting systems are shown in Fig. 4(d)-(f). Deep in the localized phase at IF = 20 
the conductivity at these frequencies is dominated by short distance single-particle hops, so the interacting and 
noninteracting systems look similar and neither show finite-size effects. At IF = 4, on the other hand, we see clear 
differences. The effective exponent a for the noninteracting system is well above one, so now the finite size effects are 
quite apparent for periodic boundary conditions. 

These considerations can be sharpened by comparing the conductivity exponents extracted from the many-body 
(MB) interacting and single-particle (SP) noninteracting data, as shown in Fig. 5 (left). We consider systems with 
open boundary conditions, and extract the conductivity exponent from a fixed frequency range that is much smaller 
than the scales J 2 /W , J 3 /IF 2 associated with short-distance hops. From our previous discussion, we expect that the 
limiting behavior for very small disorder is a ~ w 3 (because of finite size effects), whereas that for large disorder 
is a ~ w 2 (as finite size effects move to much lower frequencies and the Mott behavior is recovered). In the single¬ 
particle case, one expects the exponent a to cross over smoothly from a disorder-independent value slightly below 2 - 
on account of the logarithmic correction in Mott’s law - toward 3, and thus to increase monotonically as the disorder 
is decreased. On the other hand, for MBL, we expect that as the disorder is decreased two competing effects occur: 
the exponent decreases towards 1 for the reasons discussed in the main text; on the other hand, it is also pulled up 
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FIG. 4: (a)-(c) Frequency-dependence of a.c. conductivity in the interacting system with periodic and open boundary conditions 
(solid and dashed lines respectively) for three values of disorder W, corresponding to the thermal phase (left), near-critical 
regime (center), and deeply localized phase (right). The lowest-frequency finite-size behavior in the thermal regime goes as 
a ~ uj for periodic boundary conditions and a ~ w 3 for open boundary conditions, as discussed in the text, (d)-(f) Frequency 
dependence of a.c. conductivity in noninteracting systems, for disorder values corresponding to those in panels (a)-(c). Again, 
solid lines represent periodic boundary conditions and dashed lines represent open boundary conditions. 


Open boundary conditions 



Periodic boundary conditions 



FIG. 5: Conductivity exponent for interacting spin chains (solid red lines) and the corresponding noninteracting chains (dashed 
blue lines) as a function of disorder W, for various system sizes. The left panel shows data for open boundary conditions 
(OBCs); the right panel, for periodic boundary conditions (PBCs). Exponents are extracted from the frequency regime 
10~ 4 <w<5x 10 -3 [OBC] and 10 -4 < u < J 2 /(41U) [PBC]. Arrows indicate the evolution of the exponent with increasing 
system size. For OBCs, the exponent crosses over from the finite-size value a = 3 to the Mott value a ~ 2 with increasing 
disorder. The crossover is nonmonotonic for interacting systems but monotonic for noninteracting systems, as discussed in the 
text. For PBCs, noninteracting systems again exhibit a monotonic crossover from the finite-size exponent a = 1 to the Mott 
exponent a = 2. However, the exponent for interacting systems drops well below the finite-size value a = 1, in the regime 
where these systems are thermal. Moreover, finite-size effects on a seem negligible throughout the MBL phase. 


toward 3 by finite-size effects. Thus, we expect it to exhibit a non-monotonic U-shaped disorder-dependence with a 
minimum near the MBL transition. These expectations are borne out by the numerical data (Fig. 5): the dip of the 
many-body exponent below 2 becomes stronger for larger system sizes, approaching the value for periodic boundary 
conditions, and thus supporting the view that the true exponent for the many-body case is disorder-dependent and 
dips well below 2. 

As a final point of comparison, Fig. 5 (right) shows the conductivity exponents for several different system sizes and 
as a function of disorder, extracted from the data with periodic boundary conditions. For these exponents, we do the 
fit over several decades of data between uj l and the microscopic scales J 2 /W. Thus, these exponents can be directly 
compared to the analogous MB ones plotted in Figs. 1 and 3. We see that for any given system size, the SP data look 
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(a) W=1 



frequency uj 




FIG. 6: Distribution of squares of dipole moment matrix elements (i.e., |(n| jSj\m}\ 1 2 * 4 ) contributing to the conductivity at 

each frequency, for disorder corresponding to thermal phase (left), critical region (center), and MBL phase (right), all with 
open boundary conditions. In (c) the arrows mark the peaks corresponding to two-site and three-site hops in the MBL phase 
(see text); we expect the n-site peak to vanish at a frequency ~ 1/W n ~ , giving rise to a “shoulder” in the conductivity. Lower 
panel: frequency-dependence of the conductivity. Note that the different regimes of behavior in a(u) can be matched with 
features in the dipole moment distribution. 


qualitatively similar to the MB data, showing a monotonic increase from the finite-size dominated exponent a ~ 1 at 
low disorder (note that a ~ 1 at the MB transition due to a completely different physical mechanism) to an exponent 
approaching 2 at larger disorders. However, unlike the MB system, the SP effective exponents are strongly finite-size 
dependent and approach a constant IP-independent value on increasing system size. On the other hand, the MB 
exponents show no system size dependence in the localized phase and converge to IP-dependent values (significantly 
less than 2 for moderate disorder), further supporting our claim that the MB exponents are not finite-size effects. 


b. Nature of transitions contributing to conductivity 

We conclude with some details about the structure of the eigenstates contributing to the conductivity. The con¬ 
ductivity always exhibits a sharp feature near u> = 1, which is due to nearest neighbor resonances. For strong disorder 
(Fig. 4c), the conductivity also develops a noticeable “shoulder” at a frequency w ~ J 2 /AW due to second-neighbor 
resonances. For frequencies above this shoulder, the dominant processes are these very short range hops and the 
power law fit does not work, not surprisingly. This can be seen by looking at the distribution of dipole moment 
matrix elements, Fig. 6, of transitions contributing to a: above the shoulder in (c), a peak appears at about unity, 
corresponding to second-neighbor resonances, which is marked by the bottom horizontal arrow. Third-neighbor reso¬ 
nances can also be seen as a feature near (3/2) 2 = 2.25, indicated by the upper horizontal arrow. Data at lower values 
of IP show more such peaks, but these “shoulders” become less and less pronounced until they disappear altogether 
for IP < 8. 
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